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Abstract 

Using the Poshl- Teller approximation, we evaluate the neutrino quasinormal 
modes (QNMs) of a Kerr-Newman-de Sitter black hole. The result shows that for 
a Kerr-Newman-de Sitter black hole, massless neutrino perturbation of large A, 
positive m and small value of n will decay slowly. 
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1 Introduction 



It is well-known that there are three stages during the evolution of the field perturba- 
tion in the black hole background: the initial outburst from the source of perturbation, 
the quasinormal oscillations and the asymptotic tails. The frequencies and damping time 
of the quasinormal oscillations called "quasinormal modes" (QNMs) are determined only 
by the black hole's parameters and independent of the initial perturbations. A great 
deal of efforts have been devoted to the black hole's QNMs for the possibility of direct 
identification of black hole existence through gravitational wave detectors in the near 
future [HI2I- The study of black hole's QNMs has a long history. Most of the studies 
immersed in an asymptotically flat space time. The discovery of the AdS/CFT [3J H] 
correspondence and the expanding universe motivated the investigation of QNMs in de 
Sitter E] and anti-de Sitter [3 El El UD] space time in the past several years. 

Most methods in evaluating the QNMs are numerical in nature. Recently, using the 
third-order WKB approximation, Cho evaluated the Dirac field QNMs of a Schwarzschild 
black hole A powerful WKB scheme was devised by Schutz and Will ^3], and was 
extended to higher orders in JH]. Konoplya Ej extended the WKB approximation to 
sixth-order and calculated the QNMs of a D-dimensional Schwarzshild black hole. Zhi- 
denko ^B] calculated low-laying QNMs of a Schwarzschild-de Sitter black hole by using 
sixth-order WKB approximation and the approximation by Poshl-Teller potential. Car- 
doso calculated QNMs of the near extremal Schwarzschild-de Sitter black hole by 
using Poshl-Teller approximation, which was proved to be exactly in the near extreme 
regime [20! ■ Yoshida [21] numerically analyzed QNMs in nearly extremal Schwarzschild- 
de Sitter spacetimes. The Kerr black hole is a more general case. It is also important to 
note that the most important QNMs are the lowest ones which have smaller imaginary 
on the astrophysical aspect and the most important spacetimes are the asymptotically 
flat and now perhaps the asymptotically de Sitter which supported by the recent ob- 
servation data. So we discuss the QNMs of a Kerr-Newman-de Sitter black hole in this 
paper. Leaver [22] developed a hybrid analytic-numerical method to calculate the QNMs 
of black holes and applied to the Kerr black hole. Seidel and Iyer [23 computed the 
low-laying QNMs of Kerr black holes for both scalar and gravitational perturbations by 
using third-order WKB approximation. Berti et.al dealt with highly damped QNMs of 
Kerr black holes in |2l^ l23j. 

In this paper, we evaluate the QNMs of Kerr-Newman-de Sitter black hole for neu- 
trino perturbation. In Section 2 we consider the massless Dirac equations for massless 
neutrino in the Kerr-Newman-de Sitter black hole and reduced it into a set of Schrodinger- 
like equations with a particular effective potential. We analyse the properties of the par- 
ticular potential in section 3 and use the Poshl-Teller potential approximation to evaluate 
the QNMs of massless neutrino in Section 4. Conclusions and discussions are presented 
in Section 5. Throughout this paper we use units in which G = c = M = 1. 
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2 MASSLESS DIRAC FIELD EQUATION IN THE 
KERR-NEWMAN-DE SITTER BLACK HOLE 

Generally speaking, neutrino is a kind of uncharged Dirac particles without rest mass 
or with tiny mass. In curved spacetime, the spinor representations of massless Dirac 
equations are [TT] 
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where P A and Q A are two two-component spinors, the operator V^g denotes the spinor 
covariant differentiation. V 
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satisfy SV^b^ci) = e AC e Bbi "where e A c an d ^bd are antisymmetric Levi-Civita symbols, 
the operator V M is covariant differentiation. 

The metric of the Kerr-Newman-de Sitter black hole in the Boyer-Lindqust coordinate 
system is 
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(2.4) 
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Here, a and Q are the angular momentum per unit mass and electric charge of the black 
hole, M is the black hole mass and A is the positive cosmological constant. 
The contravariant component of metric tensor is 
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Choose the null tetrad as follows: 
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The above null tetrad consists of null vector, i.e. 



0. 



The null vector satisfies the following pseudo-orthogonality relations 



-m^mr 1 = 1, 



0. 



They also satisfy metric conditions 

g^u = l^Ji v + n^l u - m^m u - m^m u . 



(2.8) 
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Set spinor basis (^ = a^, in which A is the spinor component index, a is the spinor 
basis index, both indices are or 1. 

The covariant differentiation V^^" 4 for an arbitrary spinor £ A can be repreaented 
as the component along the spinor basis , i.e. 



where d a ^ are ordinary spinor derivatives, PI • are spin coefficients 



(2.11) 



Now let 
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m^dn = 5, d X Q = m^dn = 5 



(2.12) 



Then the Dirac equations (j2.1|) and ()2.2|) can be rewritten as four coupled equations 

(D + e - p) F x + (5 + 7T - a) F 2 = 0, 

(A + fi - 7) F 2 + (5 + (3 - r) F x = 0, 

(D + e* - p*) G 2 - (8 + 7r* - «*) Gi = 0, 

(A + fi* - 7 *) G 2 + (5 + (3* - r*) G x = 0, (2.13) 



3 



where F 1; F 2 , Gi, G 2 are four- component spinors with Fi = P°, F 2 = P 1 , G\ = Q 1 , G 2 = 
Q°. a, P, 7, e, ir, p, r etc. are Newman-Penrose symbols, while a* , (3* etc. are, respec- 
tively, the complex conjugates of a,f3 etc.. The Newman-Penrose symbols are 
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Take three transformations as follows: 
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d = e^'e im ^i(r, 0), G 2 = e-^V m ^ 2 (r, 0), 

U 1 (r,9) = g*f 1 (r,6), U 2 (r,6) = f 2 (r,8), 
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Eqs. fl2.13|) turn into 
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By using separation of variables, Eqs. (|2.18|) become 
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where Ai, A2, A3, A4 are the separation constants, and 

Ai = A3 

Then Eqs. ()2.20|) can be written as 
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Substituted \/2A with A, V2R_i with Eqs.(j222I) can be written 
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Eqs.flTZHJ) reduce to 
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Introducing the tortoise coordinate transformation from the radial variable r to the 
tortoise coordinate r* which is given by 

d A r d 



dr* c<j 2 dr' 



where 
We set 

where 
and 

Eqs. fl2.25|) is reduced to 
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eq. (f2~3!2|) and (jZHSj) changed to be 
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From eq. ()2.35|) and ()2.36|) . we obtain the radial wave equation 
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where 
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A in Eqs. (j2.23j) - (j2.38j) is a separation constant and a is the QNMs of the black hole. 
We derived the expression of A at the limit of small cosmological constant A and slow 
rotating black hole in Appendix^! It can be written as a function of I (or j), m, a and a. 



3 PROPERTIES OF MASSLESS DIRAC FIELD EF- 
FECTIVE POTENTIAL 



From the Schrodinger-like equations in eq. ()2.37j) . we can evaluate the QNMs. The 
form of V + and VI shown in eq. ()2.38|) are super-symmetric partners derived from the 
same super-potential W [12]. Ref [T3] has proved that potentials related in this ways 
have the same spectral of QNMs. Thus we deal with eq. ()2.37)) with potential V+ in eval- 
uating the QNMs. The effective potential also depends on o. 

The QNMs are decided by the effective potential. Here we analyze the dependence 
of the effective potential on parameters m, I, Q, A, a and a. The effective potential as a 
function of r is plotted for some configurations of m, I, Q, A, a and a in FigOHl From 
these figures, we can see that the dependence of V on I, A, Q is stronger than on m, 
a and a. Because there exists cosmological a constant A, the space time possesses two 
horizons: the black hole horizon r = r e and the cosmological horizon r = r c . While r 
varies from r e to r c , the effective potential V reduces to zero. 

Because of the rotation of the Kerr black hole, the azimuthal degeneracy of magnetic 
quantum number m is destroyed. On the FigHJ we show the dependence of effective po- 
tential on m, and find that negative m will increase the maximum values of the effective 
potential and decrease the position of the peak. 

In the Figl^J we show the dependence of the effective potential on angular momentum 
number I. It is clear the peak value and position of the potential increase with I. 

FigEl shows the dependence of the effective potential on the electric charge Q of 
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Figure 1: Variation of the effective potential for massless neutrino with A = 0.01, a = 
l,a = 0.1, Q = 0.1,1 = 5, j = §,m = §, §, — §, — §- 
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Figure 2: Variation of the effective potential for massless neutrino with A = 0.01,(7 = 
l,a = 0.1, Q = 0.1, m = \,l = 1,2,3,4,5, j = 1-1/2. 
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Figure 3: Variation of the effective potential for massless neutrino with A = 0.01, a 
l,a = 0.1,m= §,/ = 5, j = §,Q = 0,0.1,0.3,0.6,0.9. 




Figure 4: Variation of the effective potential for massless neutrino with a = l,a 
0.1,m = §,Z = 5,j = \,Q = 0.1, A = 0,0.01,0.02,0.05,0.1. 
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Figure 5: Variation of the effective potential for massless neutrino with A = 0.01, a 
1,Q = 0.1,m = f,/ = 5, j = §,a = 0,0.05,0.1,0.15,0.2. 
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Figure 6: Variation of the effective potential for massless neutrino with A = 0.01, a 
0.1, Q = 0.1, m = \,l = 5,j = l,u= 1,10,20,100,200. 
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black hole. The electric charge of black hole will increase the the peak value, decrease 
the position of the peak and change the position of black hole horizon. 

FigJUshows the dependence of effective potential on the cosmological constant A. In- 
creasing of A reduces the peak value of the effective potential, decreases the cosmological 
horizon radius r c and increases the black hole horizon radius r e . 

In the Fig|SJ we show the dependence of the effective potential on a. For the positive 
value of m, rotation of the black hole reduce the peak value and increase the position of 
the peak. 

For the limit of small a and the perturbation method used in Appendix El the effec- 
tive of rotation to the separation constant A is small. When we consider the dependence 
of the effective potential on a, we neglect the change of A. For A = 5, the FigJB] suggests 
the potential varies slowly as a increases and approaches a limiting position. 



4 MASSLESS NEUTRINO QNMS OF A KERR- 
NEWMAN-DE SITTER BLACK HOLE 



In this section, we evaluate the QNMs by using Poshl- Teller potential approximation 
instead of more popular WKB approximation. Konoplya used sixth-order WKB 
approximation to calculated the QNMs of a D-dimensional Schwarzshild black hole and 
compared with the result of third-order WKB approximation. Zhidenko jlHj calculated 
QNMs of a Schwarzschild-de Sitter black hole by using sixth-order WKB approximation 
and the approximation by Poshl- Teller potential. The results of Ref [Tj\ and ^Hj show 
that for large I and A Poshl- Teller potential approximation can give results agree well 
with the sixth-order WKB approximation. 

Proposed by Ferrari and Mashhoon [T3|, the Poshl- Teller potential approximation 
method use Poshl- Teller approximate potential 



It contains two free parameters (Voand b) which are used to fit the height and the second 
derivative of the potential V(r*) at the maximum 
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The QNMs of the Poshl- Teller potential can be evaluated analytically: 




(4.3) 
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where n is the mode number and n < I for low-laying modes. From eq. (|4.3jl . the real parts 
of the QNMs are independent of n. This relates to how to approximates the effective 
potential. 

The effective potential also depends on a. This will complicate matters in eq. ()4.3|) 
because there are a dependence on both sides of the equation. Thus we cannot obtain 
a directly. For slowly rotating black hole with little cosmological constant A, we can 
expand separation constant A as series of a <C 1, and expand the effective potential in 
power series of a. 

Firstly we express the position of the peak of the effective potential as series up to 
order a 5 , 

Tmax = r + no + r 2 a 2 + r 3 a 3 + r 4 a 4 + r 5 a 5 

= r + S , (4.4) 

and 

= F'(w) = V'(r ) + Z V"(r ) + ^ZlV'"(r ) 

+iEgyw (r ) + ^(5) (ro) + _L S 5 r(6) (ro) j 

(4.5) 

where r is the position of the peak of effective potential for the non-rotating black 
hole. Eq. (j2.38|) for the non-rotating black hole case, show that the expression of effective 
potential V does not depend on a and can be solved independently. We evaluate the 
coefficients r$'s order by solving this equation. The expression of r max contains a and 
unknown a. We also expand a as a = o"o + o\ a + o%c? + a^a 3 + 04a 4 + 05a 5 and plug in the 
expansion for r max , and then expand the derivation of the potential Vq performed with 
respect to r*. We plug all these expansions back to eq. ()4.3|) and solve the coefficients Uj's 
self-consistently order by order in a. 

Now we evaluate massless neutrino QNMs of a Kerr-Newman-de Sitter black hole by 
using Poshl- Teller potential approximation for A = 0.01, Q = 0.1, and a = 0.1, plot the 
results in FigJ7|and list them in Appendix EI The results show that the real parts of the 
QNMs increase with I and the magnitude of the imaginary parts increase with n. 

Figures of QNMs varying with a is plotted in the Fig|H]and plot the real and imaginary 
parts of QNMs as a function of a on the FigEl We vary a from to 0.2 to satisfy the 
condition a < 1. Because of the spherical symmetry of a non- rotating black hole, the 
QNMs is the azimuthal degenerate which can be verified in Fig|H]and FigEl a = means 
non-rotating case, which is Reissner-Nordstrom-de Sitter black hole here. Different values 
of m have the same QNMs in the FigJH] and FigEl They also clearly display the split of 
azimuthal degeneracy as a increasing from to 0.2. 

From the FigEl we see that for a rotating black hole, the real parts of QNMs is also 
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Figure 7: QNMs of massless neutrino for A = 0.01, Q = 0.1, a = 0.1. 
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;ure 9: Real and imaginary parts of the massless neutrino QNMs as a function of a for 
= 0.01, Q = 0.1,1 = 3, j = |,m = ±|,±|,±|. 
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related to the n though we use the Poshl- Teller potential approximation. For a slowly 
rotating Kerr-Newman-de Sitter black hole with little value of cosmological constant A, 
the separation constant A can be written as eq. (|A.17|) and the function of a which relates 
to n. As to the increasing of a, the real parts of QNMs increase for positive m and 
decrease for negative m, while the split of real parts for negative m is bigger than the 
positive case with the same magnitude of m. The split of real parts for n increase with a 
and the magnitude of m. n increases the real parts for the positive m and decreases the 
imaginary parts. The real parts split of n for positive m is bigger than imaginary case. 
For the negative m, the larger magnitude of m change the imaginary parts more. The 
rotation increases the magnitude of imaginary parts and the split of different values of 
m. For the positive m case, it is more complex but the FigJS] shows that the imaginary 
parts of different positive m for different n trend to the same values and the tendency is 
more clearly for large n. 

We plot the image of QNMs varying with Q in the FigEl and plot the real and 
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Figure 10: QNMs of massless neutrino vary with Q for A = 0.01, a = 0.1, 1 = 3,j = 
§,m = ±§. 

imaginary parts of QNMs as a function of Q on the FigJTTl The real parts of QNMs 
increase with Q, the split of different values of n decrease first and increase later. For 
example, while m = — | the real parts of n = is larger than n = 2 for Q = and the 
other way round for Q = 0.9. The split of imaginary parts of different m for the same n 
increase with Q. 

In Fig^J and Fig^] we plot image of QNMs varying with cosmological constant 
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Figure 11: Real and imaginary parts of the massless neutrino QNMs as a function of Q 
for A = 0.01, a = 0.1,1 = 3, j = \,m = ±|. 
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Figure 12: QNMs of massless neutrino vary with A for Q = 0.1, a — 0.1, 1 — 3, j — |, m — 
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A and the real and imaginary parts of QNMs as a function of A. Fig0] shows that A 
influence the effective potential more than other parameters. The real parts of QNMs 
and the magnitude of imaginary parts decrease with A. Just like varying with Q the split 
of different values of n decrease first and increase later. The split of imaginary parts of 
different m for the same n increase with Q. For a sufficient large A, the imaginary parts 
split of different m will be larger than that of different n, which means that in the limit 
of the near extreme A term for a slowly rotating black hole, the imaginary parts of the 
QNMs are mostly determined by m rather than n. This is why in FigEH the lines of 
m = | and m — — § seem to approach two dots for different m. 




Figure 13: Real and imaginary parts of the massless neutrino QNMs as a function of A 
for Q = 0.1, a = 0.1,1 = 3,j = \,m= ±§. 



5 CONCLUSIONS AND DISCUSSIONS 

We have evaluated low-laying massless neutrino QNMs of a Kerr-Newman-de Sitter 
black hole by using Poshl- Teller potential approximation. We adopt a further approxi- 
mation by making perturbative expansions for all the quantities in powers of parameter 
a. 

In general, the real parts of QNMs increase with I. The magnitude of the imaginary 
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parts increase with n. A character feature of the QNMs of rotating black holes is the split 
of the azimuthal degeneracy for different values m. This is clearly displayed in Figsd |H1 
IH1 For a Kerr-Newman-de Sitter black hole, massless neutrino perturbation of large A, 
positive m and small value of n will decay slowly. 

We can expand this methods to others black holes. All these works will enrich our 
knowledge about QNMs of different kind of black holes and give direct identification to 
distinguish the kind of black holes, through gravitational wave detectors in future. 
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ANGULAR EIGENFUNCTIONS AND EIGEN- 
VALUES 



The angular equations ()2.26j) can be combined into 



Aj£tS_i 



X 2 S_ 



0. 



(A.l) 



and a similar equation for S + i obtained by replacing 9 by ir — 9. For small value of 
cosmological constant A, we expand eq. (|A.l|) as series of A, say by first-order expansion 
and 



where 



v 2 



-a 2 cos 2 9V + X> 2 A + A 



| cos 2 i 



d 2 d 

d9^ + COte d9 sin 2 9 

laojm — au cos 9 — a 2 u 2 sin 2 6, 



S_i 

2 



0. 



2o . d 2o9 n 

— a coso'smfr— am a ucos9 

3 d9 3 3 

9 4 1 

42*2/1 H 2*2/i 

— a oj sin 9 H — a urn H — a sin 9 
3 3 6 

2 1 

— a 3 com cos 2 9 H — a 4 u; 2 sin 2 9 cos 2 9 

3 3 

— a 3 cu cos 9 sin 2 9 a 2 cos 2 9. 

3 2 



(A.2) 



(A.3) 
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The operator T> has no dependence on a and the solution for S_i of equation 

V S_i=-E 2 S_i (A.4) 

2 2 

can be written in terms of the standard spin- weighted spherical harmonics (Newman and 
Penrose 1966 [2Ij; Goldberg et al. 1967 [2H]) 



S_i (cos 9) e im * = _iY jm (9, <p) , E 2 = [j 



(A.5) 



where j = ^ 2 ^ with positive integer I. In general, the function s Yj m (9,ip) are defined 
in the following ways: 



s Y jm (9, if) 



2J + 1 

47T 



DL sm (ip,9,0) 



x 



(2j + 1) (j + m)\ (j - m)\ 
4tt (j + s)\ {j-s)\ 

E 



sm 



2.7 



irmp 



3 ~ s 

n 



j + s 
n + s — m 



'-iy~ s - n cot 



2n+s—m 



(A.6) 



where D J _ sm ((p, 9, 0) are the elements of the matrix representations of the rotation group 
which satisfy 



= ^2 ( 3ihtoto\3to + to) ( 3i32m 1 m 2 \jm 1 + m 2 ) D 



j 



1+m >,/3, 7 ), (A.7) 



da ^r mi («^.7)^ m2 («,A7) 



2ji + 1 

where ( j 1,72^1 ^2 li^i + ^2) are the Clebsch-Gordon coefficients and 

3 = ji + 32, ji +32 ~ 1, • • • , - J2I, 

2tt /•2tt 











da d'y sin (3 d(5. 







The function s Yj m (9, ip) satisfies the parity and the orthogonality relations 

s Y* m {9,ip) = {-l) s+m -W^v). 



dO = / / sm 



'7T /*27r 



(A.8) 

(A.9) 
(A.10) 

(All) 
(A.12) 

(A.13) 



.7^=0 



19 



With the customary definition, 



s Y jm (6, <p) 



(2j + l) (j-m)\ 



imip 



An (j + m)\ 

the operators Ct s and C +s , for au = 0, are 'raising' and 'lowering' operators: 



(A.14) 



(do Tfl CSC 9 S COt 6) s P jrn = - [(j - S) (j + S + 1)] 2 s+lPjrn 

{d e + m esc 9 + s cot 0) a P jm = + [(j + s)(j - s + 1)]* ..jP^. 



(A.15) 
(A.16) 



For small values of a and A, we can view H' = V>\ + (|a 2 cos 2 6>TAo + ^2) A as a 
perturbation operator and obtain the results by using ordinary perturbation theory, say 
by the first-order expansion 



A 



3 + 



P, ( ^a 2 cos 2 9V + V 2 ) A 



1 



-jm ) + 



(A.17) 



here 



— -jm 



a 2 cos 2 9V + V 2 A 



— -jm 



iY* m ( V l + [ -a 2 cos 2 eV + V 2 )A) _iY jm dQ 



'Idoj m - ( -r/'L - m — -a 2 m 2 | A 



a V + 1 2 V _ 1 2 ] A 
3 6 



+ 



1 2 / ■ 1 
3° l' + 2 



auo + -a 3 u;A 



2 1 

s 2 
a 00 m a 

3 2 



-jm 



A 



-jm 

sin 2 
1 



-jm 



-jm 



-jm 



cos 



-jm 







) {-\ jm 


COS 6? 



-jm 



2 2a 

'3° A 



1 4 2a / 1 ■ 

+ -a a; A ( — jm 
3 \ 2 



1 3 A / 1 ■ 

~3° 2 Jm 



cos 9 sin # (9fii 
sin 2 9 cos 2 6* 
cos sin 2 9 



■-jm 



(A.18) 



(A.19) 
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--jm 



cos 9 



cos 2 9 



.-jm ) = ( jlmO\ jm) 



J 2 



. 1 \ m 

Jo ) = 



2 ju + iy 



(A.20) 



--jm 



1 2 

- + -( j2m0\ jm) 



1 2 [3m 2 -j(j + l)] [|-J(J + 1)] 
3 3 (2j-l)j(j + l)(2j + 3) 



1 



-jm 



sin 



1 



-jm 



2 2 
2" "7 ~3~3 

2 _ 2[3m 2 -j(j + l)] [f-iO' + l)] 

3 3 
1 



( j2m0\ jm) 



j'2-0 
J 2 



1 

J 2 



1 

J 2 



(A.21) 



-jm 



(2j-l)j(j + l)(2j + 3) 
(■( >s~ /9 sin ~ —\jm 



(A.22) 



2 2 . . / 1 



35 



( j4m0| jm) 



J 2 



J; 



_2_ 2 [3m 2 -j(j + l)] [f~j(j + l)] _ 8 (2j + 1)(2j-4)! 

15 21 (2j-1)j(j + 1)(2j + 3) 35 ' (2j + 5)! 

x {(j + m)(j + m - l)[(j + m - 2)(6j - 8m + 3) - 9(j - m)(j - 3m)] 

-(j - m)(j - m - l)[9(j + m)(j + 3m) - (j - m - 2)(6j + 8m + 3)]} 



x 



3 + 



3 - 



1 



2 

2 
5 
m 



1 

~ 2 
cos sin 2 9 



9 J + 



(6j- 
J + 



9 j 



1 



-jm 



( jlmO| jm] 
1 



1 

jl-0 
J 2 



1 



2/ 5 
8(2/ - ()(2j -3)! 
5 



( j3m0| jm) ^ j3-0 



j ~2 
(6j + 7) 



(A.23) 



5j(j + l) 5 (2j + 4)! 
[(j + m) (j + m — 1) (4j — 5m + 1) — (j — m) (j — m — 1) (4j + 5m + 1)] 



3 ~2 



J "2 



, . 7 
47 + - 
J 2 



1 



-jm 



sin 9 cos 9de 



-im 



3 + 



(j2m0\ jm) ( j2 



(A.24) 
1' 



J; 



+ vK J "0( J + ^) (j2mO|jm) ( j2 ^ 1 

3m 2 -j(j + l) 21 



3 2 



±3(3 + 1) 



(A.25) 



B QNMS OF KERR-NEWMAN-DE SITTER BLACK 
HOLE 



Table 1: QNMs of massless neutrino for A = 0.01, Q = 0.1, a = 0.1 



n 



m— 7T 



m- 



m 



m- 



m- 



1 

2 





1 

1 
2 

1 
2 
3 

1 
2 
3 
4 



n 



0.1851 
0.3723- 
0.3729 
0.5561- 
0.5564 
0.5568 
0.7399 
0.7401 
0.7404 
0.7407 
0.9238- 
0.9239 
0.9241 
0.9243 
0.9246 



-0.10022 
0.094532 
-0.28322 
0.093022' 
-0.2789* 
-0.46462 
0.0925H 
-0.27752 
-0.46232 
-0.6469z 
0.092282 
-0.27682 
-0.46122 
-0.64562 
-0.82982 



0.1748- 
0.3639- 
0.3634 
0.5481- 
0.5478 
0.5475 
0.7320 
0.7318 
0.7316 
0.7314 
0.9159- 
0.9158 
0.9156 
0.9154 
0.9152 



0.099092 
0.09442i 
-0.28372 
0.0930H 
-0.27922 
-0.46562 
0.0925H 
-0.27762 
-0.46282 
-0.64822 
0.092282 
-0.27692 
-0.46162 
-0.64632 
-0.83122 



0.3792-0.094742' 
0.3808-0.28322 
0.5632-0.093102 
0.56403-0.27892 
0.5651-0.46402 
0.7471-0.092542 
0.7476-0.27742 
0.7484-0.46202 
0.7492-0.64612 
0.9310-0.092302 
0.9314-0.27682 
0.9320-0.46112 
0.9326-0.64512 
0.9333-0.82882 



0.3535-0.094542 
0.3521-0.28532 
0.5388-0.093102 
0.5379-0.27982 
0.5371-0.46732 
0.7232-0.092562 
0.7226-0.27792 
0.7219-0.46372 
0.7213-0.64992 
0.9073-0.092322' 
0.9069-0.277H 
0.9064-0.46212 
0.9058-0.64742 
0.9053-0.83292 



0.5694-0.093202 
0.5707-0.27902 
0.5726-0.46382 
0.7536-0.092602 
0.7544-0.27752 
0.7557-0.46182 
0.7572-0.64552 
0.9377-0.092332 
0.9383-0.27692 
0.9393-0.46102 
0.9404-0.64472 
0.9415-0.82812 



m— — 



m 



m= -j_ 



g - 

m =J _ 



m= — 



3 





0.5279-0.093362 














1 


0.5266-0.28112 














2 


0.5256-0.46992 












4 





0.7132-0.092692 


0.7595-0.092672 


0.7020-0.092952 










1 


0.7123-0.27852 


0.7606-0.27772 


0.7008-0.27952 










2 


0.7113-0.46512 


0.7625-0.46192 


0.6995-0.46722 










3 


0.7104-0.65232 


0.7645-0.64522 


0.6987-0.65542 








5 





0.8979-0.092402 


0.9439-0.092382 


0.8875-0.092542 


0.9497-0.092442 


Oi 


^761-0.092772 




1 


0.8973-0.27742 


0.9447-0.27702 


0.8867-0.27802 


0.9507-0.27712 


0. 


8750-0.27882 




2 


0.8964-0.46292' 


0.9460-0.46102 


0.8855-0.46412 


0.9523-0.46112 


0. 


8736-0.46582 




3 


0.8955-0.64892 


0.9476-0.64452 


0.8844-0.65092 


0.9543-0.64452 


0. 


8725-0.65352 




4 


0.8948-0.83512 


0.9492-0.82752 


0.8836-0.83792 


0.9564-0.82722 


0.8717-0.84142 
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